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Abstract 

We report results of a Monte Carlo simulation of the cj)'^ quantum 
chain. In order to enhance the efficiency of the simulation we combine 
multigrid simulation techniques with a refined discretization scheme. 
The resulting accuracy of our data allows for a significant test of an 
analytical approximation based on a variational ansatz. While the 
variational approximation is well reproduced for a large range of pa- 
rameters we find significant deviations for low temperatures and large 
couplings. 



1 Introduction 



The physics of one-dimensional quantum systems has attracted considerable 

attention, both experimental and theoretical, for a long time. Among the 
methods to treat these systems analytically the variational approximation^"^ 
has been shown to be a very powerful and useful one.^ Since, however, it 
is hard to give precise inherent error estimates for the variational approach 
it is therefore desirable to check the method against independently obtained 
data. For the 0^ chain which has been investigated both as a classicaF"^ and 
as a quantum system^'^° apparently no such independent data exist up to 
this date. 

Precise Monte Carlo simulations of many-particle quantum systems based 
on a path-integral representation of the partition function would provide just 
such an independent approach for these systems. The difficulty here is to 
achieve sufficient accuracy. Standard path-integral simulations suffer from 
well-known draw backs, such as appreciable systematic errors due to the 
necessary discretization and severe slowing down in the continuum limit. 
In order to overcome these problems a Fourier Monte Carlo simulation was 
tried some time ago for the closely related sine-Gordon chain. Even though 
preliminary data seemed to reproduce the variational approximation^ quite 
well similar results for the 0^ chain were not obtained. Unfortunately, a full 
account of these investigations was never published. A disadvantage of the 
method used in Ref.^^ is that it is not based on importance sampling which is 
a problem particularly for unbounded potentials such as the 0^ double well. 

In view of these difficulties it is therefore gratifying that recently some al- 
gorithmic improvements developed for spin systems and lattice field theories 
could successfully be transferred to path-integral simulations.^^ Multigrid 
simulation techniques'^ in particular have been shown to eliminate slowing 
down in the continuum limit for one-particle systems.'^ It seemed therefore 
worthwhile to investigate whether these algorithmic improvements may now 
render path-integral simulations of quantum chains sufficiently accurate to 
allow for a significant comparison with the variational approximation. In 
this letter we will report simulation data for the 0^ quantum chain obtained 
by combining multigrid simulation techniques with a refined discretization 
scheme. It will be shown that the accuracy of the data does allow for a 
qualified judgement about the validity of the variational approximation. 
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2 The model and variational approximation 

The system we are going to discuss is defined by the partition function 

N 

with a Hamiltonian 
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where the potential is given by 



TV 



(3) 



Here /3 = I/UbT denotes the inverse temperature, 0^ = d(pi/du, and V[{4>i{u)}] 
is the usual path-integral measure. The partition function describes a set of 
N harmonically coupled oscillators of mass Aa separated by a distance a, 
with each oscillator moving in a double-well potential. As usual, we assume 
periodic boundary conditions, 0o = 0Af- 

Following the notation of Ref.^ we introduce dimensionless parameters 
and define a coupling constant Q — hui/Eg which controls the quantum 
character of the system by determining whether the kinks are "heavy" enough 
to be treated semiclassically. In our simulations we will fix the energy scale 
by setting the energy of the classical static kink Eg = {2/3)AaujQUJi = 1. We 
also introduce the parameter R = ujq/uji which measures the length of the 
classical kink in units of the lattice spacing a. The reduced temperature will 
be denoted hj t = kBT/Eg. 

The variational approach for one-dimensional quantum systcms^'^'"'' starts 
from a quadratic trial Hamiltonian. The parameters in this trial Hamiltonian 
are determined by optimizing the Jensen-Pcicrls inequality for the free en- 
ergy. A numerical solution of the resulting set of A^(A^ -|- l)/2 self-consistent 
equations is extremely comphcated. Therefore only the hmiting cases of high 
and low temperatures and for small coupling Q have been treated in the lit- 
erature. For the latter case, which seems to be the most useful one, the final 
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result reads^ 
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where i = A(l-3L>),c:'f 



;i-3L'),andL> = E 



(cthF^-l/F^; 



with Ffc = ^hjjJkjl and o;^ = 4a;o sm^(/c7r/iV) + (J\. 

In the sequel the thermodynamic observables of interest will be the inter- 
nal energy per site u = U /N = (^F — and the specific heat per site c 

given by c = C/N = jj^- More precisely, we will be interested only in the 
anharmonic contribution to these quantities. For the free energy this is given 
by dF = F — -Pharmon = F — {1 / P) J2 ln(2sh F^). Therefore the last two terms 
of eq.(4) give the anharmonic contribution to the quantum free energy after 
subtracting the corresponding classical contribution Fdass = ln(2Ffc). 
In order to obtain analytical data for a comparison with our Monte Carlo 
results we therefore have to evaluate the configurational integral (4). 

This can be achieved by employing standard transfer integral techniques/ 
Here we have to find the eigenvalues of the transfer integral equation associ- 
ated with eq. (4) . In the thermodynamic limit — > oo only the lowest eigen- 
value survives but in order to control finite-size effects we need to compute all 
eigenvalues. With decreasing temperature more and more eigenvalues have 
to be taken into account. In particular, we observe that for low temperatures 
the two lowest eigenvalues are almost degenerate. 



3 Simulation techniques 



The partition function (l)-(3) was discretizcd Tising the Takahashi-Imada 
scheme. The discretized version of the partition function here reads 
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(5) 



where the potential is given by 

L n2f.2 N / OTA \ 2 

Here k denotes the additional index for the Trotter discretization at each site. 
The convergence of this discretization is of the order where e = hfi/L and L 
is the Trotter number. The standard, quadratically convergent discretization 
scheme is recovered by dropping the second term in eq.(6). 

Since for local update algorithms we expect a quadratic slowing down in 
the continuum limit of large Trotter numbers L^^ we applied a multigrid W- 
cycle with piecewise constant interpolation-*^^ at each site along the Trotter 
direction. Note that since we are approximating the continuum limit only for 
the Trotter discretization we do not need to apply two-dimensional multigrid 
coarsening. Also note that the interactions between the oscillators enter the 
multigrid coarsening only as constants for the polynomial expression for the 
energy on the finest grid. 

The observables of interest are the internal energy and the specific heat. 
As to the question of energy estimators we remark that the discretized par- 
tition function (5) gives rise to a so-called kinetic estimator of the energy 
Uk'^ by differentiating (U) — —dlnZ/dp ^ [4 where C/k denotes the simple 
arithmetic mean over N^n measurements of C/k in the Monte Carlo process. 
Applying a simple scaling argument one can find a different but equivalent 
energy estimator f/^,^^ based on the virial theorem with different variance. 
In order to reduce the variance of the energy estimation we may then use a 
linear combination of these two estimators. In doing so it should be noted 
that the optimal combination of the two estimators has to take into ac- 
count the individual variances and the covariance of the (blocked) individual 
estimators. Note that the energy estimators differ for the standard dis- 
cretization scheme and the improved discretization since the correction term 
in Vti is /^-dependent. For the evaluation of the anharmonic contributions 
the discretization error was further reduced by subtracting the exact values 
for the harmonic contribution at finite Trotter number L. For the standard 
discretization this improvement was already made use of in Ref..^^ Since 
we are only dealing with Gaussian integrals the exact harmonic contribution 
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can, however, also readily be found for the Takahashi-Imada scheme. ^'^ A 
full account of the simulational details discussing various systematic algo- 
rithmic refinements of path-integral Monte Carlo simulations will be given 
elsewhere. ^° 

4 Results 

We have performed simulations of the partition function (l)-(3) using the im- 
proved discretization scheme (5), (6) for different values of N and t. The 
parameter R was kept fixed aX, R — h for all simulations. The update was 
performed using a multigrid W-cycle with piecewise constant interpolation in 
Trotter direction at each site with single-hit Metropolis updating and ni = 1 
pre-, 77-2 = postsweeps. For each data point we have measured the internal 
energy using the optimally combined estimator with A'^^ = 200000 measure- 
ments taken every second sweep, i.e. rig = 2, after discarding Ue x 1000 
sweeps for thermalization. The Metropohs acceptance rates were adjusted 
to be ~ 40 — 60% on the finest grid and the same step width was used for 
all multigrid levels. The specific heat was measured by simple numerical 
differentiation of the "combined" estimator which was reweighted in a tem- 
perature interval of dt — 0.0001. These estimates gave consistent values with 
direct measurements of the specific heat using the estimators obtained by 
differentiating the discrete partition function but (slightly) smaller errors. 
All statistical errors were computed by jackkniving the data on the basis of 
500 blocks. 

Comparing the jackknife error with the canonical variance of the individ- 
ual measurements we find that the integrated autocorrelation time for both 
the kinetic energy estimator and the virial estimator never exceeded a value of 
r™*/ne < 2. Within these bounds we noticed that the autocorrelation times 
tended to be larger for low temperatures and large coupling constants. This 
observation is also reflected in the fact that the acceptance rates were roughly 
constant on all levels for high temperatures and small couplings but tended 
to decrease for lower t and larger Q. We conclude that in our simulations the 
measurements of the energies were more or less statistically decorrelated. 

Figure 1 shows the measured anharmonic contributions to the internal 
energy per site for Q — 0.1 and t — 0.1,0.2,0.3, and 0.4 as a function of 
the number of oscillators A^. Here the Trotter number L was set to L = 16 
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Figure 1: Finite-size dependence for the measured anharmonic contributions 
to the internal energy per site. Sohd hues show the variational approximation 
for finite A'^ and dotted lines show the thermodynamic limit. 



for all temperatures. The solid lines show the variational approximation 
for finite A*", and the dotted horizontal lines show the corresponding values 
in the thermodynamic limit oo. We see that the Monte Carlo data 

fully confirm the variational approximation within the statistical uncertainty. 
For high temperatures the finite-size effects arc quite appreciable but die off 
rapidly with increasing N. For low temperatures on the other hand the finite- 
size data approach the thermodynamic limit rather slowly but the absolute 
values differ only by a small amount from the asymptotic value. 

Let us now look at the temperature dependence of the internal energy. 
Figure 2 shows the measured anharmonic contributions to the internal energy 
per site as a function of the temperature t for various couplings Q. Here the 
number of oscillators was N — 300 except for t — 0.05, 0.30, 0.35, and 0.40 
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Figure 2: Anharmonic contributions to the internal energy per site as a 
function of the temperature t for various couphng parameters Q. Sohd hues 
show the variational approximation for N — 300. Dotted lines show the 
variational approximation ior N — oo. 

where we simulated a chain of = 200 oscillators. The Trotter number was 
L = 16 for t > 0.20, L = 32 for t = 0.15, L = 64 for t = 0.10, and L = 128 for 
t = 0.05. Regarding a comparison with the variational data we observe that 
the approximation again is fully confirmed for high temperatures t and small 
couphngs Q. For lower t we still find a satisfactory agreement if we also take 
into account finite-size corrections. The situation is different, however, for 
low temperatures and large couplings as can be clearly seen in Fig. 2. Here 
we find significant deviations from the variational approximation. Note that 
the error bars for the data are smaller than the data symbols. Let us take a 
closer look at the lowest temperature which we have investigated, t — 0.05. 
For Q — 0.1 and Q — 0.2 the variational approximation is confirmed within 



7 



our statistical error estimate. But already for Q = 0.3 we find a statistically 
significant discrepancy between our measured value of Rdu = 0.01641(21) 
and the variational approximation which predicts a value of Rdu — —0.01512 
for both N = 200 and N — oo. This discrepancy increases if we go to larger 
couplings. For the worst case, Q = 0.6, the variational approximation gives 
a value of Rd.u = -0.05898 for iV = oo and Rdu = -0.05982 for = 200. 
The simulation on the other hand yields a value of Rdu = —0.07447(26), i.e., 
the variational approximation deviates from the Monte Carlo results by 56 
statistical error bars even if finite-size corrections are fully taken into account. 

In order to check whether for the Monte Carlo data systematic errors 
due to the discretization would still play a role we have checked our data 
for t = 0.1 against simulations with smaller Trotter numbers L = 16 and 
L — 32. For Q — 0.1 and Q — 0.2 we found no significant differences but 
we did observe finite e effects for larger couplings Q. Their size, however, 
was small enough and in view of the fact that our discretization converges 
with the fourth order in e we believe that the remaining discretization error 
for small t and large Q is at most of the same order as the statistical errors. 
In any case, we observed that going to larger Trotter number would push 
the values down, i.e. would further increase the difference to the variational 
approximation. 

We conclude that our data differ significantly from the variational ap- 
proximation for large Q and small t. The question then arises whether these 
discrepancies are a consequence of the low coupling expansion or rather in- 
herent to the variational approach at this level. In view of the fact that 
the data fit quite well even for large Q at high temperatures it seems more 
likely that the latter is the case. On the basis of the validity of the Wigner 
expansion Giachetti et al.^ suggested that their expansion be valid as long as 

t > — ln(8i?) ^ 0.1468g^ (7) 
Stt 

This means, for Q = 0.3 we have to compare t = 0.05 with 0.013212 to 
explain a discrepancy of almost 5 statistical error bars. For our worst case of 
Q — 0.6 eq.(7) reads explicitly t ^ 0.053 and the discrepancy of 56 statistical 
error bars for t — 0.05 is indeed due to a violation of this condition. Looking 
at Fig. 2 for the largest coupling, Q = 0.6, we conclude that condition (7) is 
violated for almost all temperatures displayed in Fig. 2 only if we take the 
to mean: "larger by more than one order of magnitude." 
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Figure 3: Anharmonic contributions to the specific heat per site as a function 
of the temperature t for various couphng parameters Q. Sohd fines sfiow tfie 
variational approximation for N — 300. Dotted fines show the variational 
approximation for = oo. 



Let us finally take a look at the specific heat. Figure 3 shows the measured 
anharmonic contributions to the specific heat per site. Again the solid lines 
arc the variational approximation for A^ = 300 and the dotted lines show the 
corresponding thermodynamic limit. Due to the fact that the estimation of 
the specific heat involves a difference of statistically fiuctuating variables the 
resulting statistical accuracy is greatly reduced compared to the estimation of 
energies. Therefore our data for the specific heat do not allow for a significant 
falsifying test of the variational approximation. For the more interesting case 
of large couplings we also see that the statistical uncertainty unfortunately is 
even increasing, in particular for low temperatures. One therefore would have 
to conclude that neither a brute force increase of the statistics appears to be 
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a reasonable way of getting more accurate data. With these restrictions we 
nevertheless do see, however, that the general trend of the quantum effects 
as computed by the variational approximation is confirmed. 

5 Discussion 

Employing refined path-integral Monte Carlo techniques we have been able 
to considerably reduce the systematic and statistical errors of a quantum 
Monte Carlo simulation of the 0"^ chain. The resulting accuracy now allows 
for a significant test of the variational approximation. For small couplings we 
find that the variational quantum corrections to the thermodynamic quanti- 
ties are fully confirmed and only for large couplings and low temperatures do 
we observe significant deviations from the exact Monte Carlo data. The dis- 
crepancies may be due to the additional approximation of the small coupling 
expansion which was used to evaluate the effective classical potential of the 
variational approximation. It would therefore be interesting to see whether 
the Monte Carlo data might be reproduced by taking into account higher- 
order corrections in the coupling parameter Q. The discrepancies increase 
both for large couplings and for low temperatures. Since it is known that 
the variational approximation works better at high temperatures the other 
possible reason for the deviations of the analytical data may be an inherent 
failure of the variational approximation itself (at this level of accuracy). If 
this should be the case it would be interesting to see whether by calculating 
the higher-order corrections to the variational approach^^ one would be able 
to account for the remaining discrepancies. 
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